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ABSTRACT 

We report analysis of an 8 hr observation of PSR B0943+10 at 325 MHz performed 
at the Giant Metrewave Radio Telescope (GMRT) in India. B0943+10 is well known 
for displaying regular sub-pulse drifting and two emission modes. We investigate the 
modal behavior of B0943+10. By reconstructing an entire 'B' mode from two con- 
secutive 'B' modes, we estimate that the pulsar spends roughly 7.5 hrs in the 'B' 
mode and about 2.2 hrs in the 'Q' mode, on average. Although the pulsar can switch 
modes within one pulse, the sub-pulse drift rate changes with a characteristic time of 
1.2 hrs over the course of a 'B' mode. Under the subbeam carousel model we find the 
drift-rate changes are produced by a 10% increase in the average number of subbeams 
and a 16% increase in the carousel circulation time. We speculate that under the par- 
tially screened gap model the increase in circulation time should be related to a small 
increase in the neutron star surface temperature. 

Key words: MHD — plasmas — pulsars: general, individual (B0943-I-10) — radia- 
tion mechanism: nonthermal — polarisation — mode-changing phenomenon 



I. INTRODUCTION 

Pulsar B0943+10 is well known as a paradigm example 
of moding and subpulse drifting in pulsars. In the 'B'urst 
mode, the pulsar displays strong, regular subpulse drift- 
ing, while in the 'Q'uiescent mode the drifting disappears 
and only hints of a mo d ulat ion feature have been found 
[Rankin fc Sulevmanoval (|2006il ]. During the 'Q' mode, 
B0943-I-10 grows dimmer and a precursor component (here- 
after: PC) appears whose peak precedes the peak of the 
main pulse (hereaft er: MP) by some 50°, as discovered in 
iBackus et all (|2010l ). 

The papers in this seriej^ have investigated at length 
the dynamic behaviors of B0943-I-10. Letting Pi be the 
pulsar period, a regular modulation feature appearing at 
~ 0.47/Pi (hereafter: fz,obs) was found. In Paper I it was 
argued that this was a first order alias of the true modula- 
tion frequency lying at ~ 0.53/Pi (hereafter: /a). Paper IV 
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argued that there are gradual changes in /a over the course 
of a 'B'-mode interval with the remarkably long time scale 
of 1.2 hrs. 

This series has investigated whether B09434-10's highly 
regular drift modulation might be understood on the basis of 
the subbeam carousel model (hereafter: CM). It was demon- 
strated that the subpulse pattern of B0943-I-10 is produced 
by a carousel of 20 subbeams rotating with a circulation 
time P3 of some 37 Pi . The projection of the carousel on a 
plane perpendicular to the magnetic axis is treated as quite 
circular. 

Physically, the carousel is thought to be produced by 
'spark'-induced columns of relativistic primary plasma di- 
rected into the open polar flux tube. Under E x B drift 
the plasma columns preces s around the magnetic ax is . This 
model was introduced by iRuderman fc Sutherland! (Il975h 
and h as been further developed since [see iCil fc SendvM 

In this paper we present a nearly 8-hr long observation 
of B0943+10 carried out using the Giant Metrewave Radio 
Telescope (GMRT) in India. This long observation has al- 
lowed us to confirm the slow changes of the pulsar's behavior 
in its 'B' mode as discussed in Papers IV &i V. 

§11 provides a detailed description of the long GMRT 
observation as well as of the older Arecibo Observatory (AO) 
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observations used below. In §111 we investigate the modal 
evolution of B0943+10 and confirm the gradual variation of 
/a during the 'B' mode with a time scale of 1.2 hrs. This re- 
sult allows confirmation of the slow 'B'-mode profile-shape 
changes with the same time scale {Paper IV & V). Overall 
we found that the 'B'- and 'Q'-mode intervals had typical 
durations of some 7.5 and 2.2 hrs, respectively. §IV inves- 
tigates the observed change in the primary modulation fre- 
quency /s. We resolve a decrease in the subpulse spacing 
(P2) of about 10% over the course of a 'B' mode. We find 
that the carousel angular velocity (hereafter: cjc [~ 2tt/P3]) 
decreases by nearly 20%, corresponding to a nearly 20% in- 
crease in P3. Finally §V provides a discussion of these results. 



II. OBSERVATIONS 

1. Observing setup 

The 325-MHz pulse sequences (hereafter: PS) were ac- 
quired using the Giant Meterwave Radio Telescope (GMRT) 
near Pune, India. The GMRT is a mu lti-element aperture- 
synthesis telescope (jSwarup et al.l[l99ll ^ consisting of 30 an- 
tennas distributed over a 25-km diameter area which can 
be configured as a single dish in both coherent and incoher- 
ent array modes of operation. The observations discussed 
here used the coherent (or more commonly call ed "phased 
array") mode (|Gupta et al.ll2000l : ISirothiall2000l ') in the up- 
per of the two 16-MIIz sidebands. At either frequency right- 
and left-hand circularly polarized complex voltages arrive 
at the sampler from each antenna. The voltage signals were 
subsequently sampled at the Nyquist rate and processed 
through a digital receiver system consisting of a correlator, 
the GMRT array combiner (GAG), and a pulsar back-end. 
In the GAG the signals selected by the user are added in 
phase and fed to the pulsar backend. The pulsar back-end 
computes the power levels, which were then recorded at a 
sampling interval of 0.512 msec. 

2. Data reduction, interference mitigation 

In order to put the GMRT observation of B0943-I-10 into 
serviceable shape, there were several interference and obser- 
vational artifacts which had to be addressed. Due to instru- 
mental effects, the baseline varied on timescales of 10 s or 
more. To cope with this, a running mean subtraction was 
performed. The on-pulse regions were gated and replaced 
by random samples from nearby off-pulse regions. 2 s and 
then 5 s running means were performed accordingly and 
subtracted from the time series. 

Electrical interference at 50 Hz was strong. To mitigate 
this, after performing the running mean subtraction, an FFT 
of about 1 million samples was performed on the dedispersed 
time series, beginning at the first sample. In the frequency 
domain, the real and imaginary components corresponding 
to the regions around the fundamental frequency and its 
harmonics were replaced by randomly selected nearby sam- 
ples. An inverse FFT was performed and the same process 
was repeated on the remaining samples in the time series. 
This process could not remove all the 50-Hz interference, due 
to its unstable frequency and spectral leakage. To account 
for spectral leakage, large regions around the fundamental 
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Figure 1. Average intensity profiles for both modes. The 'B' 
mode (in black) is much brighter than the 'Q' mode (in red). 
The 'B' mode profile is taken from several thousand pulses in 
the first 'B' mode apparition (Bl). There is a ^ Scr detection of 
the PC lying ~50° before the MP during the 'Q' mode, agreeing 
with iBackus et al.l 1I2OIOI '). While the MP is dimmer during the 
'Q' mode, it is significantly wider. 



frequency and its harmonics would have needed to be re- 
moved, but this would have interfered with the pulsar signal 
and distorted the time series. 

Finally, an apparent gradual shift in the phase of the 
pulsar of about 10° during the 8-hr observation was present 
in the data. 2500-pulse averages were performed and the 
center of the MP profile was determined for each section of 
the data. The period used to bin the data into individual 
pulses was adjusted over the course of the observation to 
account for the small changes in phase of the profile. 

3. Overview of reduced observation 

B0943-I-10 is not very bright at 325 MHz, and so the in- 
creased sensitivity available at AO made the observations 
reported in Papers I & IV significantly higher quality than 
the long observation presented here. This prevents certain 
analyses performed in previous papers, but the great length 
of this observation opens a new path for inquiry. 

We confirm the presence of a precursor (PC) lying 52° 
before the main pulse (MP) observable in the total intensity 
only during the 'Q' mode (see Fig. [T}. After averaging sev- 
eral thousand 'Q' mode pulses, the PC is visible at a peak 
level of about 5 times the off-pulse rms. As with all previ- 
ous observations, the MP is brighter in the 'B' mode, and 
dimmer but wider in the 'Q' mode. 

The most important aspect of this observation is its 
sheer length. The observation spans almost 8 hr, with 23,160 
pulses available (after truncation due to running mean sub- 
traction). There were four breaks in the observation. After 
30-min the telescope was stopped for about 13 min in or- 
der to recheck the phasing. Over the following 7 hr there 
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were three ~10 min breaks in the observation spaced more 
or less evenly during which time the telescope was pointed 
to a bright nearby source to rephase the antennae. Figure [2] 
shows the pulse sequence for this long observation as a stack 



of 100-pulse averages. Figure 2(a) shows the 30 min of ob- 
servation before the long 13-min break. Figure |2(b)| shows 
the remaining ~7 hr following the break. 

The length of the modes in B0943-I-10 is on the order of 
2 or more hours, so it is highly unlikely that an entire mode 
would be missed due to such breaks, although the bound- 
ary between two modes could occur during an observational 
break. Luckily, as can be seen in figure |2(b)[ the two ob- 
served mode changes did not occur during breaks, allowing 
a good estimation of the position of the modal boundary. 

Because of the importance of the long 7-hr stretch af- 
ter the initial 13-min break, we will follow the convention of 
setting pulse to the beginning of this stretch. The modal 
sequence in this observation is as follows: 'B' mode from 
pulses -2409 to 10270; 'Q' mode from pulses 10270 to 15250; 
then back to 'B' mode from pulses 15250 to 23209. The ob- 
servational breaks fall at the following pulse intervals; —788 
to 0; 3,375 to 4,250; 10,705 to 11,100; 17,775 to 18,175. 



III. GRADUAL MODAL EVOLUTION 

Paper IV announced the discovery of behaviors associated 
with the evolution of the modes in B0943-I-10. At 'B' mode 
onset, gradual changes in the relative amplitudes of the lead- 
ing (Al) and trailing (A2) components of the MP were de- 
tected, with the ratio of the amplitudes A(2/l) decaying 
exponentially from 1.75 to about 0.2 (at 327 MHz). The 
carousel circulation time (hereafter: P3) at 'B'-mode onset 
was found to increase by about 5% as an exponential of the 
form [1 — exp{—t/T)], where r is the characterstic time. We 
argue for a larger increase of some 16%. For both the change 
in A(2/l) and P3, the characteristic time was found to be 
some 73 minutes, or 1.2 hours, (see Figs 5 and 8, Paper IV) 
Paper V found that the fractional linear polarization of the 
MP increases from 3 to 50% over the course of a 'B' mode. 

However, a full justification for these results remained 
elusive. Such dynamics can have important implications for 
understanding the origins of modes in pulsars, as long as it 
can be shown that these changes occur consistently and can 
therefore be linked with the evolution of a mode. Paper IV 
only had access to observations about 2 hours long at the 
longest, which introduced several obstacles: an entire mode 
was never observed from start to finish; only one 'B' to 'Q' 
mode transition had been seen (and that observation was 
only 986 pulses long); and the longest the pulsar had been 
observed from 'B' mode onset onward was some 3900 pulses 
(1.2 hr). Paper V only had access to short observations 4 to 
15 min in length and determined the dependence of the frac- 
tional linear polarization as a function of time after 'B' mode 
onset from the apparent time dependence of P3 established 
in Paper IV. 

The decrease of the amplitude ratio A (2/1) and the in- 
crease of P3 throughout the 'B' mode both depended on an 
imprecise argument involving the splicing together of data 
from several observations by noting the apparent exponen- 
tial nature of these changes at 'B' mode onset. 

We confirm the evolution of the frequency modulation 



during the course of 'B' mode as determined in Paper IV. 
To make the timescale of the changes determined in this pa- 
per independent of aliasing and model, we examine the fre- 
quency of the observed primary modulation feature (/3,06s), 
rather than P3 as in Paper I-V. Under the carousel model, 
P3 = N/fs^true where N is the number of subbeams and 
.f'i.true is the un-aliased f requency of the prim ary modula- 
tion feature. Paper I and I Backus et all (|201Q ) argued that 
the observed modulation is a first order alias of the true 
feature such that fa^tme = 1 — fa.obs, and that N = 20. Ac- 
cordingly, P3 = 20/(1 — /3,oi,s), so the timescales of the small 
changes in P3 will be reflected in the timescales observed in 
f3,obs given the gradual exponential behavior of the changes. 

Figure 3(a) shows a plot of fs^ots vs pulse number for 
the long January, 2010 GMRT observation, fs^obs was de- 
termined by performing 512-point Longitudinally Resolved 
Fluctuation Spectra (hereafter: LRF, see Paper I) every 50 
pulses. LRF's with a poorly defined primary modulation fea- 
ture were discarded. Each LRF overlaps the next by 462 
pulses, so adjacent points are largely redundant, but allow 
the stability of the feature to be determined. 

We confirm the characterstic time. We find that r — 
73 ± 6 min (~1.2 hr) for the increasing /3,06s, which is con- 
sistent with the findings of Paper IV. Fitting an exponential 
function to the circulation time as function of pulse number 
(see Fig. [31 for the second 'B'-mode apparition, we find 



,/3,o6s(i) ~ 0.471 - 0.022e"'/" ± 0.002/Pi 



(1) 



where t = at 'B'-mode onset. While we find that P3 
changes more than was reported in Paper IV (see §IV), 
eq.(Il]) actually agrees with /3,06s as measured in Paper IV. 
Significantly, B0943-I-10 was observed for 7960 pulses after 
'B'-mode onset, more than twice as long as in any previous 
observation, and the exponential behavior continues for the 
entire time. 

It should be stressed that using an exponential of the 
form of eq.(IT]) to describe /3,06s is not physically or theo- 
retically motivated. There is no model which predicts such 
a functional form. The primary purpose of such a fit is to 
describe the timescale of the changes in /3,06s and the am- 
plitude of the relative changes in /3,06s in a reproducible 
fashion. 

With the ~4 hours of observation before 'Q'-mode onset 
we see very clearly that at the end of a 'B' mode, /3,06s 
has increased to about 5% higher than at 'B' mode onset. 
/3,06s is largely constant during the first 'B' mode apparation 
which is to be expected at the end of 'B' mode if /3,06s 
follows the form of eq.(IT}. It is also observed that /3,06s 
has decreased 5% following the 'Q' mode, implying that the 
relation in eq.(IT]) is indeed a regular feature of the 'B' mode 
and is not an illusion due to observations spaced over the 
course of days or years. 

We are now prepared to sketch out the long-term evo- 
lution of the modes. From this long observation it is clear 
that /3,06s changes exponentially as a function of the form of 
eq.© during a single 'B' mode. If we consider this change to 
be continuous, then we may reconstruct an entire 'B' mode 



from the data within this long GMRT observation. Fig. |3(a) 



shows /3,06s as a function of time as observed. Fig. |3(b) 
shows /3,06s as a function of time for such a reconstruction. 
To reconstruct a 'B' mode from start to finish as in 



Fig. 3(b) the begining of the second 'B' mode apparition 
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(a) Pulse seqeunce before break (1621 pulses) (b) Pulse sequence after break (23209 pulses) 

Figure 2. Pulse sequences for the January, 2010 GMRT observation at 325 MHz. (a) shows the first 1621 pulses of the observation, (b) 
shows the remaining 23209 pulses which were observed following a 13-minute break in the observation. The horizontal magenta lines in 
the center panel of (b) mark the beginning and end of the 'Q' mode. The mode changes were found by examining pulse shape changes 
and the strong modulation features which are only present in the 'B' mode. Due to scintillation and de-phasing of the instrument there 
are apparent changes in brightness which do not coincide with the actual modal boundaries. The white lines mark the three occasions 
where the telescope went off-source to rephase the antennae. The center panel shows a stack of 100-pulse averages, with time increasing 
along the -|-y-axis. The left panels show the integrated flux for each 100-pulse average and the bottom panel shows the profile averaged 
over all pulses in the window. 27° around the MP are displayed. The intensity scale and the 0° longitude are arbitrary. 



(B2) was set to t = 0. /s.obs at the beginning of the first 
'B' mode apparition (Bl) is greater than at the end of B2, 
indicating that Bl comes later in the modal evolution than 
does B2. Thus, the first 'B' mode (Bl) was shifted in time 
to come after B2. 

To determine a reasonable spacing between the end of 
B2 and the beginning of Bl in the reconstructed sequence, 
an exponential function of the form of eq.([T]) was fitted to 
the /a, 063 in both regions for spacings between and 55,000 
pulses (enough to be essentially infinite). The fit is mostly 
constrained by the more steeply varying values in B2, so 
the constants in eq.([l| are not affected significantly by the 
spacing (A). The rms of the residuals was found to decay 
exponentially as rms = a -\- he'^^^" where the character- 
istic spacing Ac = 3250 pulses. Thus, the quality of the fit 
remains nearly constant for A ^ 4000 pulses. 

This sets a lower limit on the length of the 'B' mode. 
We find the 'B' mode should be at least 24,600 pulses (7.5 
hr). We can review t he plausibility of this as a lower limit, 
as well as a rough estimate of the actual 'B' mode length 
by considering the fraction of time B0943-f 10 is observed to 
be in each mode. We have 31.2 hrs of observations available 
composed of 24.0 hr in 'B' mode and 7.2 hr in 'Q' mode. 
The pulsar is therefore observed to be in the 'B' mode 77% 
of the time and in the 'Q' mode the remaining 23%. 

If the average 'B' mode length is 7.5 hrs, then we can ex- 
pect the average 'Q'-mode length to be 2.2 hrs, which seems 
to agree with observation. In the GMRT observation, the 
'Q' mode is 1.5 h rs long; whereas in th e only other complete 
'Q' mode known (jRankin fc Stapper3l2008i '). the duration is 



almost exactly 1 hr. In no other observation is an entire 'Q' 
mode seen, but there are two other instances where more 
than 1.5 hrs of continuous 'Q' mode is observed. That no 
other complete 'Q' modes are observed is unsurprising if 'Q' 
modes average 2.2 hr, given the ~2.5 hr that B0943+10 is 
visible to AO. 

It should be stressed that these timescales are estimates 
and are not expected to be exact, yet they do provide in- 
sight. Significant deviation from these values for a given 
mode should be expected. 

We may now re-address the findings of Papers IV & 
V. We confirm that the exponential changes in fa^obs as a 
function of time persist throughout the 'B' mode, thus al- 
lowing the estimation of the location of given observation 
in the 'B' mode evolution by measuring its /s.obs (or P3 as 
in Paper IV and V). Determining the exponential decay in 
the ratio of the trailing and leading components of the pro- 
file [A(2/l)] was dependent on knowing that fs^obs indeed 
evolves exponentially with the same timescale found in this 
paper. Finally, the increase in fractional linear polarization 
as a function of location in 'B' mode depended on estimating 
the position of short observations in 'B' mode from fs^obs- 

While we do not have polarimetry for this observation 
and cannot directly confirm the previously reported changes 
in polarization, we can investigate the profile shape changes. 
Figure |4] shows a stack of average profiles for the recon- 
structed 'B' mode [Fig. 3(b) along with the ratio of the 



amplitudes of the trailing component to the leading compo- 
nent (A2/A1) as a function of pulse number. Due to strong 
RFI and the lower sensitivity of GMRT compared to AO, we 
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(a) Observed primary modulation frequency (f^ oj^) for the 
2010 GMRT observation of B0943+10. The pulsar is in the 
'Q' mode approximately from pulse 10270 to pulse 15250 — 
indicated by the vertical blue stripe and labeled Ql — where 
no regular modulation is detected. In the first 'B'-mode ap- 
parition (Bl), fs^obs fluctuates around 0.47 cycles/Pi. At the 
onset of the second 'B' mode (B2), ^bs increases from '~0.45 
cycles/Pi to ~0.47 cycles/Pi. Using an exponential fit (eq[TJ, 
the characteristic time of this increase is 73 min, or 1.2 hr. 
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(b) h.oba for the 2010 GMRT observation of B0943+10 as in (a), 
but with the sections rearranged. Pulse is set to the beginning 
of region B2, corresponding to pulse 15250 in (a), and the region 
Bl is placed to begin at about pulse 11,000. An exponential was 
fitted to the data (see eq[TJ and is plotted as a solid black line. 
The 7.5 hrs from 'B' mode onset at the beginning of B2 until 
'Q' mode onset at the end of Bl represents a good lower limit on 
the length of 'B' modes. The fit appears to slightly underestimate 
fs obs S't the beginning of Bl and overestimate it at the end of 
Bl. This may indicate slight deviation from exponential behavior 
in f'i^obs- However, as is indicated in the caption below, adjacent 
points carry largely redundant information. A strong modulation 
feature at around pulse 12000 which lies around 0.475/Pi may be 
skewing adjacent points to be farther above the fit line. 

Figure 3. Observed primary modulation feature (/3 as a function of time for the Jan 2010 GMRT observation, (a) The entire 8-hr 
observation, (b) Entire reconstructed 'B' mode from the data in (a), ^bs wa-s calculated as follows. 512-point LRF's were performed 
every 50 pulses. Adjacent LRF's are largely redundant, but can help indicate the stability of a modulation feature. A sine function was 
interpolated around the bin with maximum power in the LRF to find the frequency of the modulation and the error on that measurement. 
LRF's where the modulation was poorly defined were discarded. Duo to significantly higher RFI during the first half of the observation, 
the error bars in Bl are larger than those in B2. 



are unable to measure A2/A1 as accurately as in Paper IV, 
however we confirm that the trailing component is relatively 
brighter at 'B' mode onset and grows relatively weaker as 
time progresses. 

By fitting a function of the form a + 6exp(— t/r) to 
A2/A1 we find a characteristic time r = 1.7 ± l.Ohrs which 
agrees within the error to the time reported in Paper IV. 
The fit is not well constrained, and there is significant error 
on our measurements, however we can confirm the general 
picture presented by Paper IV that the leading component 
gradually dominates the profile with a timescale on the order 
of one hour. 

We can therefore describe the evolution of the modes 
in B0943-I-10. At 'B' mode onset, the average profile shape 
chang es abruptly, a precu rsor component virtually disap- 
pears [Backus et all ||2010')], and a strong primary modu- 
lation feature becomes visible in LRF's. Over the span of 
some 7.5 hr the frequency of the feature changes as an ex- 
ponential function of the form of eq.(Il| by some 5%. The 
profile shape changes exponentially with the same time scale 
as f3,ob3 {Paper IV). The fractional linear polarization in- 
creases throughout the 'B' mode {Paper V). At 'B'-mode 
cessation the profile switches abruptly and the strong mod- 
ulation feature disappears. The PC brightens to detectable 
levels and the pulsar then spends around 2.2 hr in the 'Q' 
mode. Remarkably this behavior is a confirmation of the 



assorted observations reported in iBackus et all (|2010l ) and 
Papers I-V. 



IV. CHANGE IN MODULATION UNDER 
CAROUSEL MODEL 

We proceed to interpret the evolution of the primary modu- 
lation under the CM . The analysis presented in Paper I and 
iBackus et J] (|2010| ) showed that the 'B' mode of B0943-I-10 
displays a rotating-subbeam carousel comprised. For three 
PS's of several hundred pulses, spanning the entire 'B'-mode, 
the number of subbeams was measured to be 20. fs.obs is ar- 
gued to be a first order alias of the true primary modulation 
frequency /s.true such that /s.true = 1 — /3,06s. We can there- 
fore find the true P3 and the carousel circulation time (P3 
here, often referred to as P4): 



P3 = 



f'i,ob 



20P3 



20 



fx 



(2) 



(3) 



In this case, P3 and P3 will evolve similarly to /s.obs: they 
will increase as an exponential function of the form of eq.([l|, 
changing by about 5% over the 'B' mode. Any fractional 
changes in P3 will be reflected in similar changes in /3,06s- 
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Figure 4. Left panel: Normalized average profiles versus pulse 
number for the reconstructed 'B' mode as in Fig. |3(a)] Notice that 
the trailing component is prominent at pulse and decays relative 
to the leading component as time progresses. Right panel: Ratio 
of the amplitudes of the trailing and leading components in the 
profile (A2/A1). The fit (solid black line) is an exponential of 
the form a + 6exp(— t/r) where the characteristic time r = 1.7it 
l.Ohrs agrees with the characteristic time found in Paper IV. The 
profiles are taken from 2000-pulse averages. Pulse corresponds 
to 'B' mode start. A2 and Al were taken as the power in the bin 
corresponding to the peak of each component. Due to the shape 
of the profile, A2/A1 is sensitive to the selection of the bin for 
A2. Errors were therefore calculated by shifting the bin for A2 by 
1 in each direction and recalculating A2/A1. 



We may now investigate what changes in the carousel 
could account for the ~5% change in P3. For B0943+10, P3 
results from drifting subpulses. Under the CM, the subbeam 
carousel rotates about the magnetic axis, causing subbeams 
to appear at different longitudes for different pulses. This 
effect causes the observed subpulse drifting. 

P3 is therefore a measure of the time between adjacent 
subbeams crossing the same longitude in the profile. Two 
variations can then account for a change in P3: a change in 
the angular separation between subbeams (as measured in 
the plane perpendicular to the magnetic axis) or a change 
in the angular velocity of the carousel (hereafter: ujc with a 
corresponding change in P3. In either case, small changes in 
the viewing geometry will not affect P3. 

It is not simple to resolve which of these effects causes 
the observed change in P3, but to understand the physics of 
moding and pulsar emission we must answer this question. 
We consider what observational evidence we have concerning 
this matter. 



1. Change in carousel angular velocity ujc 

If ujc decreases, then the subbeams will cross a given lon- 
gitude on the polar cap less frequently, thus changing the 
observed P3. Similarly, the rate of change of the position on 
the polar cap of a given subbeam will decrease. This will 
cause the subbeam to appear to drift less quickly across the 
pulse window. 




i 



ID 12 14 16 18 20 

Longitude (deg) 



Figure 5. Sample results of a modfold: the track of a subpulse 
across the pulse window for the MJD 53491 AO observation, per- 
formed over pulses 1-256. Longitude in deg. is plotted on the 
horizontal axis and time in units of P\ is plotted on the vertical 
axis. A 256-pulse LRF was performed to determine P'j,^oha- The 
same pulse sequence was folded at Pz^oha a-nd plotted twice. A 
bounding box (red and black) was drawn around the subpulse 
track. Gaussians (magenta) were fitted to each row in the bound- 
ing box and a line (magenta) was then fitted to the peaks of the 
Gaussians, from which m (see eq. |4]l was determined. 



The drift rate of B0943+10 is fast enough that an in- 
dividual subpulse cannot often be observed in consecutive 
pulses: after being observed in one pulse, the subpulse has 
left the pulse window by the next pulse. To observe the track 
of a subpulse across the pulse window, we follow the method 
used in Paper I. An LRF is performed on a 256-pulse seg- 
ment of the pulse sequence from which the frequency of the 
primary modulation feature is determined. The 256 pulses 
are then folded according to the phase position within Ps.ots ■ 
This results in a single sequence of length Pa^obs from which 
we can resolve the rate at which a subbeam crosses the pulse 
window. To mitigate wrapping effects, the resulting sequence 
is plotted twice. This process shall be referred to as a mod- 
fold. Fig. O shows the results of a modfold. 

If uJc alone determines the variation in P3, P3 must 
change by some 5%, and there will be a corresponding 
change in the slope of the subpulse traverse (see Fig. [SJ- 
Under the CM, we expect the traverse of a subbeam to de- 
viate from a line due to our curved sightline (as projected on 
the polar cap) and the curvature of the carousel. In the case 
of B0943-I-10, our sightline traverse is tangential: essentially 
cutting the edge of the cone, so across the ~20° pulse win- 
dow we expect any deviation from linearity to be minimal, 
as observed. 

We may therefore treat the subpulse traverse as linear. 
Since we are concerned only with changes in the subpulse 
traverse, and since the shape of the traverse is determined by 
our viewing geometetry which is assumed to be constant, the 
linear approximation is further justified. As such, the slope 
(m) is then determined by Ps.ots and the spacing between 
adjacent subpulses in the pulse window (P2): 

Ps.obs 



P2 



(4) 
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Figure 6. Drift-band slope in periods/" longitude as a function of 
time since 'B'-mode onset for eight AO observations — as in Fig. [5] 
The position within each 'B'-mode observation was determined 
by measuring fs,obs and fitting to eq.((lj. There appears to be a 
gradual increase in the slope (m) of about 5%. 



For a constant P2, a 5% change in P^^obs will produce a 5% 
change in m. 

Fig. [5] demonstrates the determination of m for one 256- 
pulse subsequence. A modfold is performed. A bounding box 
is drawn around the track of a subpulse. The subpulses are 
quite Gaussian in form, so a Gaussian is fitted to each row 
of intensities within the bounding box. We then have a stack 
of Gaussians whose peaks align with the centers of the sub- 
pulses. The peak positions of the Gaussians are fitted to a 
line, and the slope of this line is taken to be m. The slopes 
so determined are fairly stable with respect to changes in 
the bounding box. 

A 5% change in m is very small given the accuracy of 
the method described above. The GMRT observation is not 
of sufficiently high quality to resolve m clearly, so the older 
AO observations were used to determine m as a function of 
time. For the observations used below (see Papers I, IV and 
V), m was determined for the first 256-pulse subsequence 
in the 'B' mode, then it was determined for a 256-pulse 
subsequence beginning 128 pulses later. Then the process 
was reiterated until the end of the mode (or observation). 

Fig- El shows m as a function of time. The m values ex- 
hibit a great deal of scatter compared to the expected 5% 
change we are trying to detect. Errors in the measurement 
of m would be expected to be of comparable magnitude: the 
extension of the subpulses is large compared to the total 
drift, their shapes deviate significantly from a Gaussian, the 
signal-to-noise ratio even after folding 256 pulses is still low 
enough to affect the fits, and the changing P3 at the begin- 
ning of the 'B' mode smears the apparent subpulse track. 

P3,obs decreases by about 5% during the course of a 
'B' mode. If changes in ujc were solely responsible for the 
changes in P^^obs, then P2 would remain constant, and m 
would vary proportionally to P^^obs- However, a 5% decrease 
in m is not observed — in fact, m appears to increase by 



about 5% (Fig. [6|. Thus, a carousel velocity change alone 
cannot account for the observations. 



2. Change in subbeam spacing 

We now explore the question of possible subbeam-spacing 
changes within the carousel. If the number of subbeams N is 
constant, then any such P2 changes will occur on timescales 
small with respect to P3. Such changes would correspond- 
ingly affect P3, but if m (eqU)) and cJc remain constant, then 
P3 will remain constant on average. P3 ~ 37Pi , so in a given 
256-pulse subsequence the carousel completes about 7 rota- 
tions. Thus P3 determined by an LRF performed over the 
subsequence will not be significantly affected. 

A change in P3 must therefore be attributed to a change 
in the number of subbeams: this will affect the average an- 
gular spacing of the subbeams. According to the current 
model of B0943-I-10, the number of subbeams is ~20 {Pa- 
per I). The smallest change in the average subbeam spacing 
must therefore be 5%, and in turn P3 will change by 5%. 
This conflicts what is observed: the primary modulation fre- 
quency appears to change smoothly by 5% over the course 
of a 'B'-mode apparition (see Fig. [3]) — not in discrete jumps. 

However, if the number of subbeams changes on 
timescales less than 256 pulses, smooth changes in the ob- 
served P3 could be explained. If the carousel has multiple 
configurations each with a different number of sparks, then 
the fraction of time it spends in each configuration will de- 
termine the observed P3. Gradual changes in the fractional 
time spent in each configuration would then cause gradual 
changes in P3. 

In order to distinguish variations in the subbeam spac- 
ing we can determine P2. As with the determination of m, 
a P2 change of around 5% could explain the observed P3 
variations, but achieving such precision is difficult. As dis- 
cussed in Paper I, P2 is around 10° of longitude. This is 
large with respect to the profile width, so it is difficult to 
observe single pulses with two subpulses: few pulses exhibit 
more than one and, in the few that do so, they fall on the 
profile wings where they are characteristically faint. Thus, 
we cannot measure P2 directly. 

We therefore follow the analysis presented in Paper I. A 
Fourier transform was performed on an unfolded time series. 
The frequency of the pulsar and its harmonics are strongly 
visible in the frequency domain. Additionally, components 
due to P3 falling at (A^ + fa) /Pi are visible, where A'^ is the 
harmonic number of the pulsar rotation frequency. The am- 
plitudes of these components peak around Nmax ~ 35 (see 
Figs [7] and [8]) . Nmax then represents the harmonic of Pi as- 
sociated with the primary subpulse-modulation frequency — 
or, in analogy with eq.Q, f^^^bs = l/(A^maxP2). P2 (given 
its units in degrees) can then be calculated as 360°/A'^maa;, 
where Nmax is the harmonic number corresponding to the 
peak amplitude. 

To determine /3, the fiuctuation spectrum is folded at 
the pulsar frequency. Fig. [7| shows an harmonically resolved 
fluction spectrum (HRF). The folded frequency domain is 
plotted in the bottom panel with the pulsar harmonics (now 
the DC component) removed. The left panel shows the am- 
plitude of the pulsar harmonics as a function of harmonic 
number. The center panel shows the amplitude of the com- 
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Figure 7. Harmonically-Resolved Fluctuation spectrum (HRF) 
performed on the first 256 pulses of the MJD 53492 AO observa- 
tion. The frequency spectrum folded at is plotted in the bottom 
panel, excluding the bins corresponding to the pulsar frequency 
and its harmonics. The intensity of the pulsar harmonics as a func- 
tion of harmonic number is plotted in the left panel. The center 
panel shows the intensity of the features lying at (Af + f)/Pi as 
a function of pulsar harmonic number N and frequency / v^here 
0</< ^. 
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Figure 8. A plot of modulation intensity vs. harmonic number for 
primary modulation feature, using the HRF computed in Fig. [7] 



ponents of the folded Fourier transform as a function of the 
pulsar harmonic number. 

The primary modulation peak in the folded Fourier do- 
main determines /a. Fig. [8] shows the amplitudes of the 
{N + f'i)/Pi components for the HRF plotted in Fig. [T] 
A cubic spline fit was performed to find the position of the 
peak. Generally, the amplitudes change smoothly, so that 
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Figure 9. P2 as a function of time in the 'B' mode for nine AO 
observations. P2 was measured by performing an HRF as in Fig. [7] 
Using the frequency bin corresponding to the primary fluctuation 
feature, the peak of the features lying at {N + f^)/ Pi was found, 
from which P2 = 360°/A„ia^. 



the splines fit the amplitudes well. There is typically a clear 
maximum. 

In this manner, P2 was determined as a function of po- 
sition in 'B' mode for the available AO observations. The 
position of each AO observation in 'B' mode was found by 
exploiting the changes in /3,06s- For each observation, /3,06s 
was found as a function of time by performing 256-pulse 
LRFs. The results were then fitted to eq.([T} and are plotted 
in Fig. 1^1 P2 appears to decrease in a manner similar to the 
increase in P3, changing by about 10% over the course of 
the 'B' mode. 

The average pulse profile changes shape with the same 
time scale as the changes in P3 and with a manner very 
similar to the observed changes in P2. We were concerned 
that the changing profile shape or the changing P3 could 
affect our measurement of P2 , thereby producing the illusion 
of a change in P2. To resolve this question, a modeled pulse 
sequences was produced to test our method of finding P2 
(see Fig. [To)) . 

The subpulses were modeled as a positively valued sine 
wave (i.e., 1 + sin(j>). The period of the sine wave was set 
to be 10° longitude corresponding to the observed P2. The 
pulses were modeled by multiplying the sine wave by a 'B' 
mode average profile taken from an AO observation with 
the off-pulse window set to 0. The phase of the subpulse 
sine wave was increased by 2ttPi/P3 from pulse to pulse to 
simulate subpulse drifting. 

To check if changes in the average profile affect mea- 
surements of P2 two methods were followed. First, 256-pulse 
average profiles were calculated from an AO observation be- 
ginning at 'B' mode onset. The pulse window over which 
the profile was calculated was incremented by one pulse 
and the process was completed for the remainder of the 'B' 
mode. In order to mimic the pulse shape changes, the pulses 
were modeled using these average profiles according to the 
method above. For the second method, a 256-pulse profile 
was calculated for the beginning and the end of a 'B' mode. 
An exponential function was used to interpolate the changes 
in each bin of the profile for the several thousand pulses in 
the 'B' mode, thus producing a smoothly changing average 
profile to be used to model pulses. 
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Figure 10. A simulated pulse sequence for B0943+10. A 50- 
pulse stack is plotted in the center panel. The average of those 
50 pulses is plotted in the bottom panel, and the flux for each 
pulse is plotted in the left panel. Subpulses were modeled by a 
positively valued sine wave. Pulses were modeled by multiplying 
the subbeam sine wave by an average profile which changed from 
pulse to pulse. For each pulse, the subpulse phase was increased 
by 27rPi/P3 to simulate subpulse drifting. 



The modeled pulse sequence was then analyzed in the 
same way as the observations to determine P2 ■ No change in 
P2 greater than a fraction of a percent was observed using 
the modeled data, from which we conclude that the gradual 
changes in profile shape cannot account for the observed 
change in P2. 

To see if the change in P3 affects measurements of P2 
we used both a constant and a changing 'B'-mode average 
profile while altering P3 as a function of pulse number ac- 
cording to eq.([T]) given the relation in eq.((2]). No change 
was detected in the measurement of P2 in the modeled se- 
quence, from which we conclude that the changing P3 cannot 
account for the observed changes in P2. 

Thus, we observe a decrease in P2 by about 10% during 
the 'B' mode. This change appears to be fairly continuous. 
Ps.ots decreases by about 5% over a similar timescale0 By 
eq.© we would expect these two effects to produce a ~5% 
increase in the slope of the subpulse track, which agrees with 
the measurements of m (Fig. [6|. We conclude that there is 
an intrinsic change in the slope of the subpulse track. 

Under the CM, this change in the slope indicates a 
change in uic- Accounting for aliasing, the observed decrease 
in P3,oi,3 corresponds to an increase in the true P3. From 
eqs[l]and[2]we see there is an increase in P3 of 4.1% during 
the 'B'-mode. Due to our tangential sight-line traverse, P2 
should be inversely proportional to the number of subbeams 

^ P3 variati o ns ca n in principle "Doppler shift" P2; e.g., see 
iGupta et all l l2004lh eq. (5). However, for B0943-I-10 here these 
effects are at the 1% level. 



uJcN 



(5) 



(6) 



(7) 



Combining these two equations tells us: 

P2 

OC — 
P3 

Given a 10% decrease in P2 and a 4.1% increase in P3 during 
the 'B'-mode, eq.Q tells us there should be a 13.5% decrease 
in ujc- Since P3 oc l/iOc, there is a 16% increase in P3 during 
the 'B'-mode. 



V. DISCUSSION 

We have presented a long, 8-hr observation of B0943-I-10 us- 
ing the GMRT telescope in India. Such a long observation 
of the pulsar is important for resolving questions about the 
behavior of the long modes in B0943-I-10. We confirm the 
gradual evolution of the primary fluctuation feature which 
corresponds to subpulse drifting in the 'B' mode. From 'B' 
mode onset until 'B' mode cessation, /3,06s varies accord- 
ing to eq.([T| with a characteristic time of r = 1.2 hr. This 
supports the conclusions of earlier papers in this series that 
during the course of a 'B' mode there are changes in the 
average profile with a similar characteristic time as well as a 
gradual increase in the fractional linear polarization during 
a 'B' mode. 

As a rough estimate, we find that 'B' modes will persist 
for around 7.5 hrs in length and 'Q' modes some 2.2 hrs. 
Significant deviation from these numbers is expected (and 
is observed for the 'Q' mode). 

The previous papers in this series have argued that the 
emission from B0943-I-10 can be understood according to 
the subbeam-carousel model, wherein 20 evenly spaced sub- 
beams rotate around the magnetic axis in about 37 Pi. 

We detect a decrease in P2 of about 10% during the 
course of a 'B' mode. Under the carousel model, this indi- 
cates an increase in the number of subbeams. Given about 
20 subbeams, there should be two m ore subbeams at the 
end of 'B' mode than at the beginning. iBackus et al.l l|2010l ) 
found three occasions, lasting some several hundred pulses, 
where the number of subbeams was measurable. These three 
occasions spanned the entire 'B'-mode. Using the analysis 
presented in § IV, P2 is found to decrease as expected for 
these occasions. For all three, A*' was found to be 20, which 
seems to contradict the fact that P2 measu r emen ts indicate 
an increase in N. However, in lBackus et al.l (|2010j), the mea- 
surement of N at the beginning of the 'B'-mode has a large 
error due to low signal-t o-noise in the LRF us ed to compute 
N. Thus the findings of lBackus et all l|2010l ) are commen- 
surate with a carousel composed of 20 subbeams at the end 
of the 'B'-mode and fewer subbeams (perhaps 18) at the 
beginninglfl 



3 See iBackus et al] ll2010h figures 1, 
ments of N. 



2, and 3 for the measure- 
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For A'^ ~ 20, one would expect discrete jumps of about 
1/20 = 5% in our measurements of P2. However this is not 
observed. Rather, the changes in P2 appear continuous, al- 
beit with significant scatter. Our measurements of P2 are 
calculated over several hundred pulses, so if the carousel 
model is correct the pulsar must switch between two config- 
urations with a different number of subbeams on timescales 
much less than several hundred pulses. This indicates that 
the pulsar preferentially spends more time in one configu- 
ration at the beginning of 'B' mode and gradually spends 
more time in other configurations. 

By measuring P2, P3, and the slope of subbeams as 
they cross the pulse window — which gives a measure of 
how quickly a given subbeam passes through a degree of 
longitude — we find that over the course of a 'B' mode P3 
increases by some 16% rather than the 5% found in Paper I. 

We may now speculate on a possible physical explana- 
tion of the gradual changes in P3. According to the par- 
tially screened gap model (PSG), the E x B drift velocity 
of the subbeams around the magnetic axis should be re- 
lated to the temperature of neutron star surface at the po- 
lar cap. In the PSG mode l , the vacuum gap as proposed by 
iRuderman fc Sutherland! (|l975l ^ is actually partially filled 
with ions (or possibly electrons). Heating of the neutron star 
surface due to bombardment by relativistic particles causes 
the ejection of charge carriers into the vacuum gap. This 
shields the potential drop AV across the vacuum by some 
factor < rj < 1 such that AV = tjAVrs, where AVrs is 
the p otential drop across a purely vacuum gap ijGil et al.l 
I2OO3I ) 



The drift velocity, a nd hence ljc, sh ould be proportional 
to AV, thus P3 cx 77"^ (|Gil et all 1200^ ). We observe a 16% 
increase in P3 over a 'B'-mode apparition, so we can say: 



P3,2 _ m 

Ps,! V2 

According to lGil et all l|2003l ): 



1.16 



(8) 



ri = l-exp[C{l-Tc/Ts 



(9) 

where Ts is the surface temperature of the neutron star, Tc 
is the critical temperature at which the potential drop is 
completel y shielded (i.e. A V = 0) and C is a constant. As 
argued bv lGil et al.ll|2003l ). should be at most a few per- 
cent less than Tc, so the argument of the exponent in eq.(|9} 
is close to 0. Performing a Taylor expansion and keeping the 
terms up to first order: 



r]^ -C{l-T,/T,) 



(10) 



From eq.® and ea. (|10|l we find that the ratio of the surface 
temperature at the end of a 'B'-mode to the temperature at 
the beginning should be: 



1.16Tc 



Ts,2 ^ 



(11) 



Now, Ts < Tc. So for 0.9Tc < T^ < T^ we find that 
1 < Ts, 2/^3,1 < 1.014. This means that according to the 
PSG model we should expect to see an increase in the sur- 
face temperature of less than 1.4% over the course of a 'B'- 
mode. B0943-I-10 has been observed in the X-ray regime 
(|Zhang et al.ll2005l ). and it may be possible with long simul- 
taneous X-ray and radio observations to detect such small 



changes in the temperature. This would be strong evidence 
in favor of the PSG model, along with the carousel model. 
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